Elliptic flow of cf) mesons and strange quark collectivity at RHIC 
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Based on A Multi-Phase Transport (AMPT) model, we have studied the elliptic flow «2 of (f> 
mesons from reconstructed K + K~ decay channel at the top Relativistic Heavy Ion Collider energy 
at Brookhaven National Laboratory. The dependences of V2 on transverse momentum pr and 
collision centrality are presented and the rescattering effect of </> mesons in the hadronic phase is 
also investigated. The results show that experimental measurement of V2 for <j> mesons can retain 
the early collision information before cf> decays and that the cf> V2 value obeys the constituent quark 
number scaling which has been observed for other mesons and baryons. Our study indicates that 
the <f> V2 mostly reflects partonic level collectivity developed during the early stage of the nucleus- 
nucleus collision and the strange and light up/down quarks have developed similar angular anistropy 
properties at the hadronization. 

PACS numbers: 24.10.Cn, 24.10.Pa, 25.75.Dw 



INTRODUCTION 



Elliptic flow in heavy ion collisions is a measure of the 
azimuthal angular anisotropy of particle distribution in 
momentum space with respect to the reaction plane 
The magnitude of the eeliptic flow depends on both initial 
spatial asymmetry in non-central collisions and the sub- 
sequent collective interactions. The elliptic flow is thus 
sensitive to the properties of the dense matter formed 
during the initial stage of heavy ion collision [2, bJ U LJ 
and parton dynamics at Relativistic Heavy Ion Col- 
lider (RHIC) energies. Experimentally, elliptic flow has 
been measured as functions of collision centrality, trans- 
verse momentum, (pesudo)rapidity and particle species 
0, E H 113 EH in 197 Au + 197 Au collisions from RHIC 
at Brookhaven National Laboratory (BNL). The experi- 
mental results of charged kaon, proton and pion [j| show 
that the elliptic flow first increases with particle trans- 
verse momentum following the hydrodynamic behavior 
and then reaches a saturation in intermediate transverse 
momentum region. More importantly, a Number-of- 
Constituent-Quark (NCQ) scaling has been discovered 
for identified particle elliptic flow in the intermediate px 
region for baryons and mesons. The data of pion ellip- 
tic flow V2 are somewhat higher than the NCQ scaling, 
which can be attributed to the large contribution to the 
pion yield from secondary decays 0, 0|. The NCQ 
scaling at RHIC is an important indication for the effec- 
tive constituent quark degree of freedom at hadronization 
and the formation of hadrons through parton coalescence 
mechanism Similarly, a Number-of-Nucleon scal- 

ing for anisotropic flow of light nuclear clusters in nuclear 
collisions at Fermi energies was recently demonstrated by 



* Corresponding author. Email: ygma@sinap.ac.cn 



Ma et al. ^(|; which is interpreted as a consequence of 
nucleon coalescence mechanism. 

Strange quark dynamics is a useful probe of the dense 
matter created at RHIC. Enhanced strangeness produc- 
tion 01 has been proposed as an important signal for 
the formation of the Quark-Gluon Plasma (QGP) in nu- 
clear collisions. The dominant production of ss pairs via 
gluon-gluon interaction may lead to a strangeness (chem- 
ical and flavor) equilibration time comparable to the life- 
time of the QGP whereas the strangeness equilibration 
time in a hadronic fireball is much longer than the lifetime 
of the hadronic fireball. The subsequent hadronization of 
the QGP is then expected to result in an enhanced pro- 
duction of strange particles. In particular, it has been 
argued that with the formation of QGP the production 
of <f> mesons is enhanced. Furthermore <f> mesons could re- 
tain the information on the condition of the hot plasma 
at hadronization because <fi mesons interact weakly in 
the hadronic matter |18|. The measurement of 4> mesons 
has been of great interest in the study of collision dy- 
namics and the properties of the dense matter created at 

RHIC mmim 

We use A Multi-Phase Transport (AMPT) model to in- 
vestigate effect of parton dynamics related to c6 mesons. 
The AMPT consists of four main components |23 : initial 
conditions, partonic interactions, conversion from par- 
tonic matter into hadronic matter and hadronic interac- 
tions in collision evolution. The initial conditions, which 
include the spatial distribution of participant matter, 
minijet partons production and soft st ring excitations, 
are obtained from the HIJING model [23. Scattering 
among partons are modelled by Zhang's parton cascade 
(ZPC) [24(, which calculates two-body parton scatterings 
using cross sections from p QCD with screening masses. 
In the default AMPT model |25j partons are recombincd 
with their parent strings when they stop interacting, and 
the resulting strings fragment into hadrons according to 
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the Lund string fragmentation mode l 1261 . In the AMPT 
model with string melting scenario |27|. a quark coales- 
cence model is used to combine parton into hadrons. The 
evolution dynamics of the hadronic matter is described 
by A Relativistic Transport (ART) model [2^. Details 
of the AMPT model can be found in pi]. 

In this paper, we present a detailed study of the el- 
liptic flow of <p mesons at the top RHIC energy based 
on the AMPT model with the string melting scenario 
22]. The string melting scenario is believed to be much 
more appropriate than the default AMPT scenario since 
the energy density in these collisions is much higher 
than the critical density for the QCD phase transition. 
The AMPT model with string melting scenario has suc- 
cessfully described the elliptic flow of stable baryons 
and mesons ^2^, |27|]. It can also describe higher-order 
anisotropic flow parameters v n including the odd-n ones 
|2*flj . In this work, we focus on the <j> mesons. We follow 
the method used in experimental analysis by reconstruct- 
ing <f> mesons in the final state from fitting the invariancc 
mass distribution of all K + and K~ pairs with a Breit- 
Wigner function including the intrinsic decay width of <j> 
mesons [3(J • The parton scattering cross section is chosen 
as 10 mb. The transverse momentum and collision cen- 
trality dependences of <f> meson v-i have been studied and 
the NCQ-scaling has also been observed in the AMPT 
model. In addition, the rescattering effect on (f> meson 
i>2 has been investigated in the hadronic phase using the 
ART model where rescattering processes of <fi mesons and 
their decay daughters are included [3lJ . We note that our 
study based on the AMPT model for cj> meson elliptic flow 
is different from that of Ref . [3^ , where dynamical quark 
coalescence model has been used for studies of (j> meson 
production and elliptic flow. 



II. 



ANALYSIS METHOD 
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turn (px), azimuthal angle ($) of pairs that formed from 
all permutations of candidate K + with K~ at a given <p 
rapidity range (|y| < 1.0). The resulting m Mra distribu- 
tions consist of the <j) signal superimposed on a large back- 
ground that is predominantly combinatorial. The shape 
of the combinatorial background was calculated using the 
mixed-event method |33j| . A Breit-Wigner function plus 
a linear function can describe the rrii nv distribution well. 
The left panel of FigQ] shows the nii nv distribution for 
<j> mesons from minimum bias collisions using ~ 1.8M 
Au+Au Monte Carlo AMPT events. 

The elliptic flow vi is calculated by the STAR stan- 
dard method based on the distribution of particle 
raw yields as a function of azimuthal angle $ with re- 
spect to the event plane angle 'J. The raw yields of cf) 
mesons are extracted from fits to rtii nv distribution in 
each pt and $ — W bin. The event plane angle W is 
determined from the azimuthal distribution of charged 
tracks within a window of 0.2 < pt < 2.0 GeV/c and 
pseudo-rapidity \rj\ < 1.0, which is used as an estimate of 
the reaction plane angle [34l l3l| . To avoid autocorrela- 
tions, tracks associated with a <f> candidate are explicitly 
excluded from the event plane calculation. The right 
panel of FigQ] shows the azimuthal angular distribution 
of raw <j> yields with respect to the event plane from the 
minimum bias collisions in the 0.4 < pt < 3.0 GeV/c 
range. Dashed line is the fitting result with a function 
d(( ^ jv ^ ) = A[l + 2v 2 cos2($-^ r )], where A is a normaliza- 
tion constant. The finite resolution in the approximation 
of event plane as reaction plane smears out the azimuthal 
angular distribution and leads to a lower value in the ap- 
parent anisotropy parameters . This event plane res- 
olution is determined by dividing each event into random 
sub-events and calculating the difference in event plane 
angles between sub-events. We obtained an event plane 
resolution of 0.91, which is 20% larger than the experi- 
mental resolution reported by STAR [HI- This is due to 
the fact that the number of tracks per event used in our 
simulation is larger than that from the data. In order to 
verify our resolution correction, we have also calculated 
the V2 with respect to the real reaction plane (\& = 0) 
which is known a priori in our model calculation. From 
Fig. [3 the V2 extracted from the true reaction plane is 
in good agreement with the one extracted from the event 
plane corrected for resolution effect. Our result also illus- 
trates that the experimental elliptic flow analysis method 
can faithfully describe the magnitude of the elliptic flow 
for d> mesons. 
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FIG. 1: (Color online) Left panel: mi nv distribution from mini- 
mum bias (0-80%) Au+Au collisions at ^sjvjv = 200 GeV. Dashed 
line is a fit of Breit-Wigner function plus a linear background. Right 
panel: azimuthal angular distributions of raw <j> yields with respect 
to the event plane. Dashed line represents the fit result. The trans- 
verse momentum range of the figure is 0.4 < pj< < 3.0 GeV/c. 



The reconstruction of <j> mesons was accomplished by 
calculating the invariant mass (rrii nv ), transverse momen- 



III. 



RESULTS 



The upper panel of Fig|3] shows the elliptic flow V2 of <p 
mesons from minimum-bias AMPT 197 Au + 197 Au col- 
lisions at y/SNN = 200 GeV. Experimental data Q of 
K + +K~ and p + p are also presented for comparison. 
We note that the <fi meson V2, in compaison with t>2 of 
charged kaons and protons all from the AMPT model, 
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FIG. 2: (Color online) Comparison of elliptic flow values obtained 
from using real reaction plane (RP) and event plane angle corrected 
for resolution effect. 



satisfies the mass ordering behavior of v 2 predicted by 
the hydrodynamic model calculation |36j | in the low pt 
region. However, the AMPT calculation results of v 2 
for charged kaons and protons are about 25% larger than 
the experimental data in the pr < 1.5GeV/c region. This 
may result from the large parton scattering cross section 
(10 mb) used in this Monte Carlo AMPT calculation. 
In order to explore the intermediate pr (1.5< pr <4.0 
GeV/c) phenomenon, where the so-called NCQ-scaling 
in elliptic flow for identified particles has been observed 
at RHIC |j| and the quark coalescence or recombination 
mechanism has been used to explain the scaling 0, 0] , 
the large parton scattering cross section is needed in or- 
der to produce the magnitude of the elliptic flow in the 
intermediate pt region matching the experimental mea- 
surement [27^ . 

Our AMPT calculation indicates that the v 2 of <fi 
mesons at the intermediate pr seems to saturate and to 
follow the same behavior as that of K + +K~ . The lower 
panel of Fig. |3| shows elliptic flow v 2 normalized by the 
number of constituent quarks for charged kaons, protons 
and cj> mesons. In the intermediate pt region oipT/n q 
> 0.6 GeV/c the elliptic flow of charged kaons, protons 
and <f> meons from the AMPT calculation seems to sat- 
isfy the NCQ scaling. This result implies that u, d and s 
quarks in the initial partonic matter formed in relativistic 
heavy ion collisions develop significant collectivity with 
strength characterized by v 2 /n q . 

We have also studied the collision centrality depen- 
dence of <j> meson v 2 at several centrality intervals: 0-20%, 
20-40% and 40-80% as well as - 80%. Fig.Hshows that 
for each centrality bin, v 2 (j>t) of <fi mesons increases with 
their px- Among different centrality bins, the values of 
v 2 {pt) increase from central to semi-peripheral collisions, 
which can be understood as a result of increasing initial 
spatial eccentricity. 
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FIG. 3: (Color online) Top panel: px dependence of the i>2 for cf> 
mesons compared with K + +K~ and p + p. Data are taken from 
Ref. |3 . Dot-dashed lines are the fitted results with function ( 

f vo (n) = T-i ; — , an i , , i , - dn, where a, b, c and d are the fit 

parameters, n is the constituent-quark number ). Bottom panel: 
NCQ scaled i>2 as a function of NCQ scaled transverse momentum. 
The error bars represent statistical errors only. 



As the K + + K~ pair from the decay of a cj> meson is 
likely to undergo rescatterings in the medium during the 
hadronic evolution, this might lead to a reconstructed 
K + K~ invariant mass situated outside the original cf> 
meson mass peak. It is thus of interest to study the 
in-medium rescattering effect in details. This is carried 
out by turning on and off the ART process during the 
hadronic evolution in the AMPT model. The elliptic flow 
v 2 of 4> mesons without ART indicates the elliptic flow 
developed before the hardonic rescattering stage. In con- 
trast, the v 2 of <p mesons after ART includes all contribu- 
tions from both partonic and hadronic stages. In Fig. [SI 
the i>2 for (j> mesons before ART are directly extracted 
from the AMPT model without the ART processes. In 
that case, mesons are explicitly present and do not 
need to be reconstructed from K + + K~ decay channels. 
The v 2 after the ART processes is reconstructed from 
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FIG. 4: (Color online) vi of <f> meson as a function of pj< in the 
centralities of 0-20%, 20-40%, 40-80% and 0-80%. The error bars 
represent statistical errors only. 
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FIG. 5: (Color online) The hadronic rescattcring effect on elliptic 
flow V2 of 0-meson from AMPT model with string melting scenario. 



K + + K~ pairs and the hadronic rescatterings are mostly 
due to kaon rescatterings in the hadronic stage. Error 
bars are statistical errors only. An interesting feature in 
Fig. [S] is that the two scenarios are in good agreement 
with each other after pt > 0.4 GeV/c. In this case, the 
final hadronic rescattering effect on the <fi meson elliptic 
flow can be ignored within the errors. Our study of <fi 
V2 confirms that <p mesons can retain useful information 
from the early stage of the nuclear collisions. 

IV. SUMMARY 



In summary, we have presented a study of elliptic 
flow of <j) mesons using reconstuctcd K + K~ pairs from 
minimum-bias 197 Au + 197 Au collisions at ^Jsnn = 200 
GeV in a multi-phase transport model with string melt- 
ing scenario. The t>2 of 4> mesons seems to exhibit a 
similar behavior as other mesons. A NCQ-scaling phe- 
nomenon of elliptic flow has been observed for <j> mesons 
from the reconstruction of K + K~ pairs. The coefficient 
V2{pT/n q ) of 4> mesons represents essentially the momen- 
tum space anisotropy of constituent strange quarks that 
have arisen from the partonic collectivity developed in 
the initial stage of heavy ion collisions. The collision 
centrality dependence of elliptic flow for <f> mesons has 
also been studied. It is found that the c/> meson elliptic 
flow increases from central to semi-peripheral collisions 
as a result of increasing initial spatial eccentricity. We 
have also studied the in-medium hadronic rescattering 
effect on elliptic flow of <f> mesons. The results confirm 
that within error bars, our reconstructed <f> meson V2 can 
retain the early information before it decays. Comparing 
our predictions with the RHIC data for the elliptic flow 
of cj> mesons is expected to shed light on these issues. 
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